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Honeycomb  Setback  acceleration 

Mitigator  Projectiles 

Impact  Gun  interior  ballistics 

Pulse  shaping 
Artillery  simulation 

ZO.  AHTNACT  (CmMmm  mt  mtmm  alNt  If  m«mi'  aW  MmMIY  *r  Mm*  mmtm) 

The  Harry  Diamond  Laboratories  4-  and  7-in.  air  guns  are  used 
to  slowly  accelerate  projectiles  and  components  contained  therein 
to  speeds  up  to  900  ft/sec.  The  impacts  of  these  projectiles  are 
designed  to  produce  pulse  shapes  (projectile  decelerations)  that 
simulate  the  interior  ballistic  acceleration  environments  occurring 
in  a real  gun.  The  continuity,  momentum,  and  energy  conservation 
equations  are  applied  to  the  impact  testing  of  projectiles,  and 
calculations  obtained  therefrom  are  oomoareil  w4i-h  ovnoT-imoni-ai  air 
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gun  data  of  projectile  time-deceleration.  Data  refer  to  i^ojectile 
speeds  of  300  to  6S0  ft/sec,  accelerations  up  to  40,000  g,  with 
energies  of  14,000  and  43,000  ft-lb  (for  the  4-  and  7-in*  guns, 
respectively)  impacting  on  aluminum  honeycomb  mitigatofs'* 'h'SVxng  a 
variable  crush  area,  uniform  static  crush  strengths  of  725,  2000, 
4000,  and  8000  psi,  and  a velocity  (strain-rate)  dependent  dynamic 
crush  strength.  . 

Within  the  limits  of  the  projectile  energy,  mass,  and  avail- 
able mitigator  strengths,  it  is  shown  that  the  entire  pulse  shape 
can  be  controlled  by  shaping  the  crush  area  of  the  mitigator  and 
by  use  of  stock-item  railroad  springs.  The  mitigator  shape  con-  j 
trols  the  pulse  rise  and  steady  parts,  and  the  springs  control  the 
pulse  fall.  Horeover,  it  is  shown  that  high  transient  stress  wave 
effects  produced  by  impact  can  be  eliminated  from  the  pulse  by 
shaping  the  mitigator  and  by  the  use  of  springs,  thereby 
attaining  a very  smooth  pulse. 
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INTRODUCTION 


The  Harry  Diamond  Laboratories  4-  and  7'-in.-diameter  x 90*£t- 
length  gas  guns  are  used  to  provide  a simulation  of  artillery  interior 
ballistic  acceleration  environments  of  ordnance  projectiles*  In  this 
simulation^  the  projectile  (called  a bird)«  having  equipment  on  board  to 
be  test-evaluated/  emerges  from  the  muszle  of  the  4»  or  7-in*  gun  with 
an  energy  of  14,000  or  43,000  £t-lb,  respectively,  and  with  speeds  up  to 
900  ft/sec*  The  equipment  in  the  bird  is  mounted  so  that  the 
acceleration  force  (pulse)  on  impact  is  in  the  same  direction  as  the 
acceleration  force  (setback)  in  the  weapon.  The  pulse  is  caused  by  the 
crushing  of  an  aluminum  honeycoinb*  mitigator,  located  between  the  bird 
and  a mxnentum  exchange  mass  (KEH),  which  is  at  rest  before  impact.  The 
mitigator  can  either  be  launched  to  move  with  the  bird  (in  which  case  It 
is  attached  to  the  bird  nose)*  or  may  be  attached  to  the  HEM  (fig.  1). 
For  a nonelastic  KEM,  the  ratio  of  KEH  to  mitigator  masses  is  of  the 
order  100,  and  the  ratio  of  MEM  to  bird  masses  is  of  the  order  10.  The 
test  simulation  is  designed  so  that  the  terminal  velocity  of  the  bird  is 
approximately  zero  and  its  entire  rixxnentum  is  transferred  to  the  HEM  or 
to  the  mitigator  and  KEM.  A description  of  the  gas  gun  facility  and  its 
operation  is  given  by  H.  D.  Curchack.^ 


The  pulse  experienced  by  the  bird  is  comprised  of  essentially 
three  parts:  rise,  steady,  and  fall.  The  rise  and  steady  parts  occur 
during  the  crushing  of  the  mitigator,  and  their  characteristic  features 
are  determined  primarily  by  the  bird  mass  and  the  dynamic  crush  strength 
of  the  mitigator.  In  addition,  particularly  at  high  crush  speeds,  the 
acceleration  of  the  crushed  mitigator,  and  thereby  the  mitigator 
density,  affects  the  bird  deceleration  when  the  mitigator  is  attached  to 
the  KEH.  As  will  be  shown,  aluminum  honeycomb  mitigators  can  be 
designed  to  provide  a range  of  smooth  rise  and  steady  decelerations, 
consistent  with  the  imoact  energy  and  mass  of  the  bird.  The  fall 
characteristics  depend  on  the  bird  and  KEM  masses  and  the  elasticity  of 
the  system,  namely,  the  elasticity  of  the  mitigator  and  the  springs 
provided  in  the  MEM.  Accordingly,  the  pulse  fall  can  also  be  varied  so 
that  a wide  range  of  pulse  shapes  can  be  obtained. 


The  conservation  equations  of  mass,  momentum,  and  energy  are  given 
for  determining  the  forces  acting  on  and  the  motions  of  the  individual 
system  components  as  functions  of  time.  A comparison  of  experimental 
and  calculated  data  is  presented  to  indicate  the  validity  of  the  theory. 


^H«  D*  Curchack,  Artillery  simulator  for  Fuse  Evaluation*  Shock  and 


Vibration  Bulletin  (Dec  1970).  Also  reported  in  Harry  Diamond 
Laboratories  TR-1330  (Nov  1966). 

*Aluminum  honeycomb  is  also  commercially  available  under  such  trade 
names  as  tubecore,  spiralgrid,  etc* 
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ttegAtJivo  No.  49-166-835  1974 
Conical  aitigator  attached  to  bird  with  aasa  MEM. 


Negative  No.  49-186-833  1974 
Bird  and  double  wedge  luitigator  attached  to  spring  MEM. 


Figure  t.  Inpact  Components 


Aluminum  honeycomb  is  cellular  in  structure  and  the  crushing 
occurs  in  ceiiuiar  columns.  Experimental  data  show  that  shear  forces 
developed  between  crushed  and  uncrushed  parts  are  insignificant  in 
con^arison  with  the  static  and  dynamic  crush  forces. 


For  a czush  at  the  bird  interface,  the  crushed  material  lies  in 
the  region  bounded  by  the  crush  front  and  the  flat  nose  of  the  bird. 
Referring  to  figures  2a  snd  3a.  the  bird  has  penetrated  the  mitigator  by 
the  distance  6 and  the  shaded  region  denotes  the  crushed  material. 
Since  the  honeycomb  structure  is  uniform  across  the  section,  the  crush 
front  is  always  planar  for  a mitigator  whose  area  remains  constant 
during  the  crush.  figure  3a.  For  the  mitigator  having  the 
tricu\gular>shaped  section  shown  in  figure  2a.  because  of  the  column-like 
crush*  the  crush  front  is  not  planar  but  is  fully  determined  by  the 
vertex  angle  a.  the  bird  penetration  6*  and  the  mass  densities  of  the 
crushed  and  uncrushed  material.  To  attain  a continuous)y  soooth  pulse, 
the  mitigator  is  designed  so  that  the  instantaneous  crush  takes  place  at 
the  weakest  section  and  progresses  continuously  toward  the  stronger 
sections.  To  this  end«  the  mitigator  section  area  is  made  a monotonlc 
nondecreasing  function  of  the  bird  penetration,  and  the  instantaneous 
effective  mitigator  crush  area  is  the  instantaneous  mitigator  face  area 
with  the  bird.  E)^eriments  show  that  the  initial  crush  force  is  about 
SO  percent  larger  for  an  uncrushed  mitigator  than  one  for  which  some 
crush  has  occurred.  Thus,  experiments  show  that  the  crush  continues  only 
at  the  shaped  end  of  the  mitigator.  even  after  the  instantaneous  crush 
area  becomes  equal  to  that  of  the  unshaped  face  (fig.  2b).  However,  the 
crush  generally  occurs  at  both  ends  of  a flat-faced  mitigator. 
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For  a crush  at  the  MEM  interface  and’  penetration  only  by  the  MEM, 
the  above  description  remains  the  same,  except  that  the  mitigator  has 
been  turned  around  (so  that  its  weakest  section  faces  the  MEM)  and  moves 
with  the  bird  in  the  gun.  figure  3b. 

In  practice,  mitigators  with  an  axisymmetric  nose  cannot  be  used 
because  the  mitigator  and  bird  center-of-mass  axes  will  not  perfectly 
coincide  and  the  bird  will  tend  to  tunble  as  it  impacts.  Consequently, 
for  stability,  mitigators  shaped  with  off-axis  nose  projections  are 
used,  figures  4.  S.  and  6. 

In  the  following  presentation,  the  length  of  the  shaped  end  in  the 
direction  of  the  bird  motion  is  called  the  altitude  length  and  the 
remaining  length  in  this  direction  is  called  the  base  length. 

The  dynamic  crush  force  of  the  aluminum  honeycomb  is  significantly 
larger  than  the  static  crush  force  emd  was  determined  empirically.  For 
a given  biid  mass,  the  rise  and  steady  parts  of  the  pulse  can  be 
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a.  Schematic. 


Negative  No.  49-186-837  1974 
b.  crush  of  conical  mitigator 


Figure  2.  Shaped  mitigator  crush. 


designed  by  shaping  the  mitigator;  i.e.«,  the  cnish  area  is  a designed 
function  of  the  bird  or  MEM  penetration  and  the  bird  speed  relative  to 
the  MEM.  The  mitigator  is  always  long  enough  to  av  .'d  complete  crushing 
or  “bottoming,”  since  the  forces  at  “bottoming"  greatly  exceed  the 
dynamic  crush  force. 


The  mitigator  retains  elastic  energy  when  the  bird  speed  relative 
to  the  MEM  is  zero.  The  amount  of  the  associated  compression 
displacement  is  small  and  thereby  the  fall  time  is  short,  typically 
about  100  usee  when  a pure  mass  MEM  is  used.  The  mass  MEM  is 
approximately  12  to  IS  times  larger  than  the  total  mass  accelerated  in 
the  gun.  The  MEM  mass  i;  selected  so  that*  together  with  the  honeycomb 
elasticity,  the  terminal  velocity  of  the  bird  will  be  approximately 


zero. 
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A spring  HEM  has  been  deviseU  tc  obtain  longer  fall  times  and  to 
allow  for  shaping  the  pulse  fall  (fig.  7).  This  KEH  consists  of  a 
forward  section,  which  uaJces  the  impact  interface  with  the  mitigator, 
two  stock'ltem  railroad  springs,  an  aft  section,  and  a washer,  bushing, 
and  connecting  bolt  to  hold  the  device  together.  The  device  is  free  to 
compress,  but  cannot  ’'fly*  apart  following  impact  because  of  the  bolt 
action.  In  a typical  shot,  the  two  railroad  springs,  each  weighing  5.3 
lbs,  provide  a total  stored  energy  corresponding  to  linear  spring 
loading  of  86,400  lbs  and  maximum  deflection  of  0.331  in.  Typically, 
this  stored  energy  provides  for  about  a 300  nsec  fall  time  vAen  the  bird 
mass  is  2000  grams  and  the  peak  bird  deceleration  exceeds  15,000  g. 
Smaller  bird  decelerations  may  not  cause  the  full  deflection  of  these 
springs  and  a reduced  stored  energy  is  available.  Other  fall  times  can 
be  obtained  by  varying  the  bird,  mitigator  and  MEM  masses,  varying  the 
sltigator  strength,  and  by  changing  the  energy  capacity  of  the  springs. 
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The  above  spring  MEM  is  about  one-halt  to  one-third  the  wjight  of 
the  pure  mass  KBH.  Accordingly,  the  bird  speed  relative  to  the  spring 
HBM  beccnes  zero  at  a larger  value  of  the  bird  speed  relative  to  the 
gun.  In  practice,  the  fall  time  la  limited  by  the  fact  that  the  final 
bird  speed  relative  to  toe  gun  auat  be  zero  or  slightly  positive  (i.e. , 
the  bird  should  not  move  baclcwtirds  toward  toe  gun  and  damage  the 
piyload).  The  acceleration-time  function  for  the  fall  is  fully 
determined  by  the  toyeical  constants. 

3.  CONSERVATION  EQUATIONS 

3.1  Mass  Conservation 

Referring  to  the  sketch  of  the  mltigator  given  in  figure  3a,  we 
assvaie  that  the  mitigator  is  attached  to  the  HEM  and  the  crush  proceeds 
from  left  to  right.  The  original  mltigator  length  is  L,  the  bird  has 
penetrated  the  mltigator  by  the  distance  £ and  toe  crushed  material  is 
contained  in  toe  shaded  space  (i-  j) . Since  the  cross-sectional  area  is 
not  significantly  affected  by  the  crush,  mass  conservation  regulres 

p/pX  - (S-l)/S  , (1) 


toere  p and  pi  are  the  mass  densities  of  toe  original  mltigator  and  the 
crushed  mltigator  in  the  shaded  region,  and  S”i/4. 

3.2  Momentum  Conservation 

In  deriving  the  momentum  equation,  we  observe  that  toe 
centers-of-mass  for  the  crushed  and  uncrushed  parts  of  the  mltigator 
move  at  speeds  which  differ  from  01  and  02  but  depend  on  01,  02,  and  S. 
Consequently,  it  appears  necessary  to  derive  a tomentum  equation  from 
the  relation  expressing  the  center-of-gravity , Y,  of  toe  entire  moving 
system  in  terms  of  the  centers-of-gravity,  YI,  of  the  system  individual 
masses,  MI.*  This  is  given  by 


MlYl  + M2Y2  + M4Y4  + M5Y5 


*The  author  is  especially  indebted  to  H.  D.  Curchack  of  HDL  for  this 
consideration,  the  starting  equation  (2),  and  helpful  discussions  that 
resulted  in  equations  3(A)  through  5(A). 


14 


’^1'^ 

IgQ 

p 

Kf-^-' 


Case  A:  Mitigator  Adjacent  To  «EH  And  Mitigator  Crushes  At  Bird 

Face. 

The  nonentun  cf  the  systeo  is  MlOO,  and  tme  differentiation  of 
equation  (2)  gives  the  jnomentum  equation 


MlUO  = MV  ” mul  + K2U2  +H404  + MSU5  + H4(V4-V5). 
Referring  to  figure  3, 


(2A) 


Y4  = (l+«>/2,  Y5  = (L+l)/2 

' 

i = U2  + 

(01-U2)S,  i « 01,  L 

Pi 

u4 » l±i 

01  <•  02  + (01-02)5 

04  2 

2 

U5  * — 

2U2  + (U1-U2>S 

2 

2 ' 

1 i 

«4  « pAA 

F M5  « pA(L-t)«  and 

fe'-  i 

M4  « pAS(01-02)  = pAU-b). 

Substituting  the  above  quantities,  we  obtain 


(MU4  + H5U5  + M4  (V4-Y5)  = + 


M4(H-S)U1  , M5(2-S)02  , M5(S)01  ^ 
2 2 2 


pAS 


(Ul-02)  ((i+6)  - (!/»•£)  1 . 


Since  pAS  (6~L)  - -M4(S-1)-M5S,  there  results 


M4U4  + M5U5  + M4(Y4-Y5)  “ H4U1  + M5U2,  and  the  momentum 


equation  becomes 


T,:i»v;^A;s 


MlUO  = J«  • (Ml+M4)01  ♦ (M2+M5)U2 


The  nooentun  equation  (3A)  Is  exactly  that  which  would  result  by 
assuming  that  the  entire  crushed  and  uncrushed  masses  move  at  speeds  Cl 
and  U2,  respectively. 

Since  no  external  forces  act  on  the  system,  time  differentiation 
of  equation  (3h)  gives 


(M1+M4)  A1  + (M2+M5)  A2  + M4  (01-02)  • 0. 


The  forces  transmlcted  to  (HltH4)  and  to  (H2tH5)  originate  at  the 
location  I,  see  figure  3a.  The  time  rate  of  change  of  mass,  M4,  gives 
rise  to  a force  acting  on  (HltH4) , so  that 


A1  - -(P+M4  (01-02)  )/{Ml+M4)  (5A) 

«4iere  F is  the  dynamic  crush  force  of  the  mitigator.  The  value  for  F is 

generally  a function  of  (01-02) , and  experiments  using  aluminum 

honeycomb  indicate  that  F is  larger  than  the  static  crush  force,  F.  In 
addition,  F varies  linearly  with  the  crush  area  normal  to  (01-02) , and 
therefore  also  varies  with  the  depth  of  penetration  for  shaped 
mltlgators. 

The  only  force  that  con  be  transmitted  to  (H24H5)  is  F, 

Accordingly , 


A2  • F/(M2+M5) 


Case  B»  Mitigator  Haves  With  Bird  And  Mitigator  Crushes  At  MEM 
Face. 

As  in  Case  A,  we  start  with  the  time  derivative  of  equation  (2) . 
Since  the  mitigator  moves  with  the  bird,  there  results 


(M1+M3)  DO  « Mf-  MIDI  + M202  + M404  + M5D5  + M4  (Y4-Y5)  . (2B) 
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langim  that  the,  entire  nase  eystem  Is  subjected  to  a yelocity  rOOf  then 
the  HEH.  and  bird  nove  to  the  left  at  sp,eeds  U2^.  • -UO  •¥  (J2  and 
Ul^‘  • -00  +’’U1.  The  ctmii  occurs  at  the  MEM  face  and  now  proceeds  from 
right  to  left.  This  crush  Is  similar  to  that  of  Case  X,  where  the  crush 
occurs  at  the  bird  face  and  the  crush  proceeds  from  left  to  right. 
Since-  V4«  YSf  M,  M4,  and  MS  are  the  same  as  in  Case  X,  we  obtain  a 
similar  result 


M4U4*  + MS05*  + M4  (Y4-yS)  - M402l>+  MSUl^  . 


HW(  superlsfioelng,  a speed  + 00  so  that  the  bird  and  mltigator  proceed 
to  the  right',  -equation  (2B)  becomes 


(M1+M3)  00  » MV-  (M1+M5)  01  + (M2+M4)02  (38) 


Thus,  the  momentum  equation  (38)  is  the  same  as  that  which  would  result 
by  assuming  that  the  entire  crushed  and  uncrushed  masses  nove  at  speeds 
02  and  01.  respectively. 

Again,  since  no  external  forces  act  on  the  system,  time 
differentiation  of  equation  (38)  gives 


(MltHS)  XI  + (H2+H4)  A2  + K4  (i)2-01)  « 0 . 


(48) 


The  force  transmitted  to  (H1+M5)  at  the  crushed  and  uncrushed 
mltigator  Interface  is  produced  solely  by  the  dynamic  crush  force  of  the 
mitlgator. 


w. 


XI  - -F/(M1+M5)  . (58) 

Accordingly,  this  force  (in  the  opposite  direction)  plus  the  force 
arising  from  the  change  in  momentum  per  unit  time  of  the  mltigator 
occurring  at  the  interface  of  the  crushed  emd  uncxushed  masses. 
M (01-02) , is  the  force  acting  on  (H2+M4) . Accordingly , 


A2  - (F+Me  (01-02)  )/(M2+M4) 


(68) 


Hi 


Cage  C;  Mitlgator  Adjacent  To  MEM  And  Hitigator  Crushes  At  HBH 
Eace«  and 

Case  Pi  Hitigator  Moves  With  Bird  And  Hitigator  Crushes  At  Bird 
Face. 


In  each  case,  stress  wave  effects  cause  undesirable,  non-smooth, 
stnnewhat  random  pulse  shape  irregularities  and  some  crush  may  occur  at 
each  end  of  the  mitigator.  Accordingly,  these  two  cases  are  not 
considered  practical  for  impact  testing. 

3.3  Energy  Conservation 

In  both  Cases  A and  B,  energy  is  dissipated  by  deformation  of  the 
mitigator.  In  each  case,  the  defomation  energy  amounts  to 


El  = /F  (01-02)  dT 


In  addition,  the  system  loses  energy  at  the  crush  front  (the  interface 
separating  the  crushed  and  uncrushed  masses)  through  collision  between 
the  crushed  and  uncrushed  mitigator  masses.  Denoting  the  mass  of  a 
particle  by  M,  the  loss  in  kinetic  energy  produced  by  this  type  of 
impact  amounts  to 


1/2  EM(U1-U2)2  , 


where  Ul  and  U2  are  the  velocities  before  and  after  intact.  ^ The 
inpacts  are  inelastic,  since  the  relative  velocity  between  the  involved 
ouisses  is  rero  followin9  inpact*  In  our  case,  the  progression  of  the 
crush  front  involves  a continuous  inpact  process  so  that  the  energy 
dissipated  by  this  mechanism  amounts  to 


1/2  /M4  (Ul-U2)2dT  . 


The  total  energy  dissipated  by  the  system  is  El+W.  Calculations 
show  that  the  ratio  H/Elso  for  UO-lOO  (where  values  of  Ul-U2<100 
ft/sec) , and  Increases  with  increasing  (U1-U2)  to  0.05  when  UO  • 500 
ft/sec. 


^E.  T.  Whittaker,  A Treatise  on  the  Analytical  Dynamics  of  Particles 
and  Rigid  Bodies,  4th  ed.,  Dover  Publications  (1944),  pp.  234-5. 
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Case  At  The  kinetic  energy  of  the  bird  prior  to  impact  is  the 
energy  of  the  system.  The  energy  balance  at  any  time  during  impact  is 
given  by 


EO  « EltE2-H)  , 


uhers  El  and  H are  given  by  equations  (7)  and  (8) , 


EO  - MJ.OoVa  and 


E2  - (M1+M4)U1V2  + <H2+M5)  02^/2  . 


Case  Bi  The  kinetic  energy  of  the  bird  plus  the  mitlgator  prior 
to  impact  is  the  energy  of  the  system.  The  energy  balance  at  any  time 
during  impact  is  again  given  by  equation  (9)  with  El  and  H given  by 
equations  (7)  and  (8) . but  where 


EO  • (M1+M3)«oV2  and 


E2  - (M1+M5)  01^/2  + <M2tH4)  02^/2  . 


4.  COKPOTER  FROGIWHS 


Programs  were  written  for  Cases  A and  B of  section  3 for  both  mass 
and  spring  hem's,  giving  rise  to  the  impacts  previously  described.  In 
each  case,  the  motions  of  the  bird,  mitigator  and  MEM  are  completely 
determined  by  the  monentum  equations  of  section  3.  The  mitigator  is 
shaped  to  produce  a smooth  pulse.  The  crush  area  of  the  shaped  face 
Increases  linearly  with  the  crush  length  from  0 to  4m  in.^  and 
thereafter  remains  constant  and  equal  to  4m  in.^>  The  dynamic  crush 
force  is  assumed  to  vary  linearly  with  (U1-U2) , and  linear  spring 
constants  are  assumed  in  the  honeycomb  and  spring  HEM. 


*The  crush  area  is  slightly  smaller  than  4m  In.^  because  the 
manufacture  of  honeycomib  requires  an  opening  of  at>out  0.5  in.  inside 
diameter. 


.^i>'jy~''^t-t'^’ ‘' ' •-  ■ ! ' ; 


P> 

■J^\ 


fc. 


!;■ 

ia-f 

I?- 


» c: 


T 


b)  tcfort  cruth 


During  crufh 


Figure  9.  Notation  for  shaped  nitigatoir  attached  to  bird  with  spring 
HEM. 


1 


f 


E ■ t 


The  fall  is  now  governed  by  a 3-spring,  4-oaoe  system  (fig.  9b).  The 
spring  constants  Zl,  Z2  axe  determined  by  input  displacements  Cl,  C2, 
and  the  forces  acting  on  Ml,  M2  at  the  time  U1-U6;  Z3  is  the  value 
appropriate  for  the  railroad  springs  in  the  MEM.  The  pulse  terminates 
when  the 'forces  acting  on  Ml,  M2,  M3,  and  M6  are  simultaneously  tero. 
Again,  Cl  and  C2  are  adjusted  to  provide  agreement  between  the 
calculated  and  experimental  pulse  fall  time. 


m 


I With  the  appropriate  equations  for  A1  and  A2,  this  program  can  be 

run  for  Case  A. 

4.3  PUI.SE3 

The  equations  and  program  describe  the  same  impacts  as  PULSEl.  In 
PULSE3,  the  experimental  pulse  fall  A1(T)  is  assumed  known  together  with 
Z1«Z2,  since  RfO.  The  C1”C2  are  adjusted  to  provide  agreeement  with  the 
experimental  pulse  fall  A1(T).  Thus,  PUI.SE1  and  PULSE3  assume  linear 
and  nonlinear  honeycomb  elasticity,  respectively.  Hence,  PULSE3  can  be 
used  to  empirically  determine  the  nonlinear  elasticity  of  the  mitigator. 

4.4  JMEM  and  JBIRD 


I 


JMEM  and  JBIRD  are  abridged  versions  of  VARYB  and  PULSEl  for  the 
mitigator  crushing  at  the  HEM  and  bird  interfaces,  respectively.  These 
programs  are  presented  to  show  that  the  momentum  and  energy  balances 
described  in  section  3 are  satisfied,  using  the  stated  expressions  for 
Al,  A2,  etc.  JBIRD  and  JMEM  do  not  take  into  account  mitigator  shaping. 
Moreover,  for  the  above  purpose,  it  was  convenient  to  neglect  the 
elasticity  terms  (which  are  point  functions)  and  the  programs  terminate 
when  U1«U2. 


?■ 
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It  is  an  easy  matter  to  include  the  JMEM  and  JBIRD  programs  within 
VABYB  or  PULSEl.  The  conservation  of  momentum  and  energy  could  then  be 
demonstrated  in  VARYB  and  PULSEl,  with  appropriate  terms  included  for 
the  springs.  However,  it  was  preferred  not  to  further  add  to  the 
printout  of  these  programs. 

5.  THEORETICAL  AND  EXPERIMENTAL  RESULTS 

The  impact  of  two  solid  bodies  results  in  the  production  of 
coiq>ression  waves  at  the  is^ct  surface.  In  a one-diaensional  analysis, 
the  waves  propagate  through  each  body  at  the  elastic  and/or  plastic 
sound  speed,  reflect  at  the  and  surface,  and  propagate  back  toward  the 
ispact  surface.  The  waves  partially  reflect  and  transmit  at  the  impact 
surface,  and  the  process  repeats  until  the  energy  contained  in  the  waves 
is  dissipated.  The  strengths  of  the  waves  are  a function  of  the  impact 
force.  For  or  steel  ispacting  bodies,  the  wave  speed  through 

each  body  is  approximately  200,000  in./sec,  so  that  the  wave  system  can 
produce  mih«t«ntlal  variations  of  the  force  acting  at  each  section  of 
each  body  in  a typical  pulse  duration  of  1 msec. 

The  data  described  in  sections  5.1  through  5.4  and  given  in  tables 
I through  III  and  figures  10  through  15  refer  to  experiments  made  in  the 
4-in.  gun  at  bird  ispact  speeds  of  400  to  500  ft/sec  and  “tigator 
static  crush  strengths  of  8000  ± 800  psi.  The  results  of  additiowl 
teats,  described  in  section  5.5  and  given  in  figures  16  ^ough  18, 
refer  to  bird  in^ct  speeds  in  the  range  390  to  650  ft/sec  and  nitigator 
static  crush  strengths  between  725  and  4000  psi* 

AS  an  example  of  the  effect  on  the  entire  pulse  caused  by  the 
svstem  of  shock  waves  arising  from  a high  initial  inpact 
the  acceleration-time  history  arising  from  the  ® 

flat-faced  aluminum  honeycomb  mitigator  moving  with  the 
and  a mass  HEM,  figure  10,  shot  nunber  1352.*  The  bird  a^ 
were  driven  by  atmospheric  air  through  a vacuum  and  attained  an  impact 
speed  of  470  ft/sec.  The  mitigator,  bird,  and  ,i:m  masses  were  °-26, 
1 38  and  24.4  kg  and  their  lengths  were  2.01  (original  uncrushed 
mitioator  length),  3.5,  and  5 in.  The  bird  was  aluminum  and  the  HEM  was 
ftee?  JirSagatir  static  crush  force  was  100,000  ± 10,000  lbs.  As 
will  be  discussed  later  in  this  section,  the  dynamic  crush  force  excee^ 
the  static  crush  force  and  appears  to  be  » ""oobooi' 

nondecreasing  function  of  the  bird  speed  relative  to  the  MEM. 


*A11  shots  refer  to  tests  made  with  the  Harry  Diamond  Laboratories'  4- 
and  7-in.  air  guns. 


22 


TABLE  U. 

StMUn  or  CALCULATED  FCtCES  AMD  MOTIONS  OF  SY8T1M  COMFOMSMT8 

fOR  CONICAL  MITXGATOR  ATZACIED  TO  BOD.  WITH  SFEIMC  MDl  AMD 

LIMEAE  BOMnCOKB  ELASTICXTT 

ntsEi 

16I83EST  11/81/79 

i n*l  Vl«;  U0«?  t2t0M«*3«460 

TIHC 

-AI 

U1  Y1  A8-  US  Y8 

A6  U6 

Y6 

r n rt 

f. 

.1 

46|.  .11  .t  .1 
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.1  .8 

• 00 

.8  .8  .8 

iie> 

7.9 

447.  .97  .1  .1 

.88 

8.8  14.8 

•0t 

38.3  4.'5  rc9 

au. 

14.6 

411.  1.17  .9  1.6 

.88 

11.9  46.6 

• 04 

99. 1 6'.'8  18'.'3 

3M. 

19.8 

394.  1.93  8.9  7.1 

.81 

9.8  81.9 

• 82 

73.'6  4'.'6  89'.'3 

4II< 

83.1 

889.'  1.98  4.9  18.3 

.88 

9.9  186.'6 

• 23 

84.'6  8.4  S9'.'3 

str.' 

84V  4 

818.  8.81  6.3  39.7 

.89 

1.1  116.9 

.37 

88.8  .7  8S>'4 

Ml. 

84.'9 

188.  8.41  6.9  8 8.9 

.13 

.7  69.6 

• 46 

89.'9  .'3  89'.'t 

674. 

84.4 

71.'  8.'9I  6.8  97.3 

.21 

1.8  71. '6 

• 52 

87.9  .8  80.7 

i SPRIMQ  CONSTANTS  CljC2>t  «0SS«*0OS 

TIME 

-Al 

Ut  XI  A2  U2  X2 

A3  U3 

X3 

A6  U6 

711. 

84.4 

49V-  .199  6.8  183. 

• 000 

138.9  186. 

.298  *-12.1  68. 

7SI. 

9.8 

83.'  .888  9.'4  118. 

.889  4-IB2.0  37.  ' 

271  a.'7  63'. 

sia.' 

13.'6 

4.  .831  4.8  188.' 

.888 

77.9  30. 

.239  -9.6'  86.  " 

851. 

>9.  .888  3.9  187. 

.'888 

.8  84. 

.194  -7.'8  48. 

911. 

.1 

•9.  .888  8.8  133. 

.888 

.0  -4.- 

.148  -5.7  39. 

991. 

.1 

•9.  .881  1.6  136. 
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.'9  -4. 
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.1 
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.0  -4. 
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A6 
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• 1 

• 1 
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5 
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Ml?  410V 
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.9 

1.6 
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1V6 

46.'9 

.84  59'. 

6 

6.1  lira 

388.’ 

19.9  354'. 

1.53 

2*5 

7V0 

• 01 

9.8 
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J.  ia  74’. 

1 
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U'92 

4.' 5 
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0 
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24V 5 206V 
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39'.'9 
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1.1 
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8 
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.3 
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4 

.8  86'.'4 

SPRINO  CONSTAHTS 
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-Al  U1 

XI 

A2 
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A3 

U3 

X3  A6 

U6 

788. 

19.7  48 

.848 
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•005 

•428 
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.384 

• 
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186'.' 
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TABLE  lA.  SUMMARY  OF  CALCULATED  FORCES  AND  HOTIMS  OF  SYSTEM  CCMP(WENT5 
FOR  CONICAL  K1T1CA70R  ATTACHED  TO  BIRD»  WITH  SPRING  MEM  AND 
LINEAR  HONEYCOMB  ELASTICITY  (Cont*d) 

PULSEI  I61UEST  ll/2t/7S 

F0»}  Vt«I  UO*?  121000iO«46e 

TIME  -At  U1  Yl  A2  U2  Y2  A6  U6  Y6  F R 
674*  24*4  71*  2*S0  fi«2  97*3  *21  1.0  71*6  >52  87*9  *0  8 


SPRING  CONSTANTS  C1«C2«7  «06S/ *001 


TIME 

-At 

Ul 

XI  A2 

U2 

X2  A3 

U3  X3 

A6 

U6 

700. 

24.4 

49* 

.065 

6*0 

103. 

•000  138.9 

ns.  *298 

♦ -12. 

1 61. 

7S0. 

10.6 

22. 

.028 

5*4 

112. 

• 001  -93.4 

38.  .271 

2.7 

63. 

800. 

13.S 

3* 

.036 

4*8 

120. 

-C00  76»6 

32.  .239 

-9.6 

57. 

850. 

• 0 

-8* 

• 000 

3«9 

127. 

•000  .0 

89.  .194 

-V.9 

43» 

900. 

• 0 

•8. 

• 000 

2*9 

133* 

•000  .0 

-4.  .142 

-S.T 

40. 

950* 

• 0 

-8. 

• 000 

1.7 

136* 

•000  .0 

•4.  .083 

-3.4 

33. 

1000. 

• 0 

-8. 

.000 

• 4 

138. 

•000  .0 

-4.  .019 

-.8 

30. 

1018. 

.0 

-8. 

• 000 

.0 

138. 

•000  .0 

-4.  *000 

.0 

30. 

TABLE  IB.  SUMMARY  OF  CALCUUTEO  FORCES  AND  MOTIONS  OF  SYSTQf  COMPONENTS 
FOR  CONICAL  MITIGATOR  ATTACHED  TO  blIU>>  WITH  SPRING  HEM  AND 
NONLINEAR  HONEYCOMB  ELASTICITY 


TIME 

-Al 

Ul 

Yl 

A2 

U2 

Y2 

A6 

U6 

Y6 

F 

R F3 

674* 

24.4 

71. 

2* 

50  6*2 

97.; 

i *21 

1.0  71.6 

.52  87.9 

.0  80.7 

Cl«  JC2»  JT3«-*  . 

• 06«  *06«  .00102 

TIME 

-Al 

Ul 

XI 

A2 

U2 

X2 

A3 

U3 

X3 

A6 

U6 

700. 

20.4 

50  • 

•057  6.0 

103 

• 060 

•6.2 

69  • 

.300 

-1.3 

71. 

750. 

12.7 

24. 

•043  5.6 

112 

• 054 

-17.9 

46< 

.277 

-3.2 

68* 

800. 

5.5 

10. 

•032  4.9 

121 

.036 

-4.4 

28  • 

.246 

-6.8 

59. 

850. 

3.8 

3 

.020  4.1 

128 

.020 

• 4 

26. 

.203 

-0.3 

48. 

900. 

3*1 

-3- 

•005  3.0 

134 

.008 

-11.0 

21. 

.151 

-3.6 

4v* 

950. 

1.9 

-6. 

.000  1.9 

138' 

.000 

.0 

6. 

.092 

-3.7 

34. 

1000. 

• 8 

-9. 

.000  .6 

139 

• 000 

.0 

6. 

.029 

-1.2 

30. 

1024* 

• 0 

-9. 

.000  .0 

140' 

• 000 

.0 

6< 

.000 

.0 

30. 

ltiS2EST  ll/2t/7S 


F0«;  UU;  U0>?  121000«O«460 
TIME  -Al  Ul  Yl  A2  U2  Y2 
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A2  U2  X2 
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•062  S«6  112* 
•048  4*9  121* 
•037  4.1  128* 

• 026  3.0  134. 

•013  1.9  137. 
•001  .6  139. 

•000  .0  140. 


TIME 

-Al 

Ul  Yl 

674. 

24*4 

71.  2. 

Cl*  ;C2*  ;T3 

■ ? .08. 

TIME 
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'SE  IB.  SUMMARY  OF  CALCULATED  FORCES  AND  MOTIONS  OF  SYSTEM  COMPONENTS 
FOR  CONICAL  MITIGATOR  ATTACHED  TO  BIRD,  WITH  SPRING  MEM  AND 
NONLINEAR  HONEYCOMB  ELASTICITY  (Cont'd) 

PULSE3  ll:4SEST  M/21/7S 

F0-1  Vl«;  U0»?  12te09««3>460 


A smooth  pulse  can  be  obtained  by  reducing  the  strength  of  the 
wave  system  arising  from  the  impact.  To  this  end,  we  first  note  that 
the  elastic  and  plastic  wave  speeds  are  at  least  an  order  of  riiagnitude 
greater  than  the  above  bird  speed  and,  because  of  energy  dissipation, 
the  strength  of  the  waves  diminishes  as  they  traverse  the  body.  Hence, 
the  entire  pulse  can  be  made  as  smooth  as  desired  by  producing,  in  time, 
a continuous  progression  of  sufficiently  small  strength  impacts,  each 
impact  increasing  the  deceleration,  until  the  desired  bird  deceleration 
IS  attained.  Physically,  this  kind  of  impact  force  can  be  impressed  on 
the  bird  by  shaping  the  mitigator  so  t-hat  the  mitigator  crush  area 
smoothly  and  continuously  increases  with  the  bird  penetration  over  the 
period  of  the  pulse  rise.  In  addition,  springs  in  the  MEM  are  helpful 
in  providing  a continuously  smooth  increasing  impact  force  with  time. 

For  smooth  pulses,  the  conservation  equations  of  section  3 may  be 
applied  to  the  determination  of  the  motions  of  the  bird,  mitigator,  and 
MCM  for  the  entire  pulse.  The  following  discussion  is  limited  to  a 

comparison  of  calculated  and  test  results  for  shaped  mitigators  with 
both  mass  and  spring  MEM's. 

5.1  Conical  Shaped  Mitigator  on  the  Bird  with  a Hass  MEM 

Typical  experimental  data  of  bird  deceleration  with  time  are  shown 
in  figures  11a  and  11b  for  shots  1703  and  1709,  using  the  conic  ^ 
mitigator  illustrated  in  figure  4.  Tae  mitigator  altitude 
section)  and  base  (constant  area  section)  lengths  were  1.79  i 
1.70  xa«,  respectively.  The  experimental  curves  are  fairly  smooth,  but 
exhibit  times  where  the  pulse  abruptly  changes  slope;  such  changes 
typically  occur  on  the  passage  of  a shock  or  stress  wave.  Thus,  in  1703 
and  1709,  stress  waves  produce  changes  in  deceleration  of  about  4,000  g. 
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Figure  10.  Setback  acceleration  for  flat-faced  mitigator  attached 
to  bird  with  nass  tCM.  Ml»1380,H2«S4400»  M3«*260^  and 
00^468.  Shot  1352. 

The  static  crush  loading  for  the  alundnum  honeyconb  in  the  above 
shots  was  e3q>eri»entally  determined  as  90»000  t 9000  lbs.  The 
experimental  data  accumulated  for  setback  deceleration  shew  that  the 
dynamic  criish  loading  F is  strain-rate  sensitive  and  is  given 
approximately  by  the  relation 

F = FO  tH-Vl(Dl/U0)l.  (12) 


Here,  A ^ 12.56  in.2  is  the  instantaneous  crush  area  (the  mitigators 
were  shaped  so  that  the  crush  area  increased  linearly  with  the  bird 
penetration  distcuice) , and  FO  and  VI  are  constants  to  be  determined  from 
the  experimental  data.  As  will  be  shewn,  the  calculations  can  usually 
be  brought  into  agreement  with  the  test  data  for  FO  approximately  S to 
20  percent  larger  than  the  static  crush  strength  and  VI  in  the  range  0,5 
to  0.3  for  (U1-U2)  < 500  ft/sec.  The  calculated  curves  were  determined 
by  using  the  con^uter  program  VAFY6  of  section  Af  (F0,V1)  pairs  were 
selected  to  provide  agreement  with  the  experimental  data. 


Specifically  f agreenent  with  esq;>eriment  was  sought  for  A1(T). 
There  is  some  elasticity  in  the  mitigator,  and  for  this  reason  the  time 
at  which  the  mitigator  crush  ends  (that  isy  when  U1bTJ2)  cannot  be 
precisely  determined  from  the  experimental  values  of  Ui  or  Al.  However, 
since  momentum  is  conserved  (eq  (3A))  and  (3B)),  the  value  of  U1»U2  is 
independent  of  PO  and  Vl«  Thus,  Ul  at  the  end  of  the  crush  is  known. 
Except  for  a possible  error  of  less  than  30  psec  in  determining  the 
beginning  of  the  pulse,  the  experimental  Ul(T)  is  a}so  well  known. 
Hence,  T at  the  end  of  the  crush  is  also  well  defined.  The  crush  time 
markedly  decreases  with  increasing  P.  Hence,  (PO,  VI)  pairs  are 
selected  in  program  VAPYB  to  provide  A1(T)  agreement  with  experiment 
consistent  with  the  above  determined  crush  time.  Although  a unique 
(FOfVl)  pair  is  not  ob'.ained,  good  agreement  with  es^eriment  is  limited 
to  (F0,V1)  pairs  over  small  ranges  of  FO  and  VI. 


The  fall  part  of  the  pulse  begins  when  Ul«U2  and  ends  when  Al»0. 
The  time  duration  and  shape  of  the  Al(T)  pulse  fall  depends  on  the 
elasticity  of  the  uncrushed  part  of  the  mitigator  and  is  independent  of 
(F0,yi). 

The  possible  strain  under  load  of  the  bird  and  MEM  is  too  small  to 
account  for  the  long  time  duration  of  the  pulse  fall.  Hence,  the  pulse 
fall  can  only  be  explained  by  elasticity  of  the  mitigator.  Accordingly, 
we  assume  equal  springs  at  each  end  of  the  mitigator  and  vary  the  values 
of  the  spring  constant  to  obtain  agreement  with  the  experiment.  The 
experimental  data  of  the  pulse  fall  are  not  sufficient  to  more  than 
grossly  define  the  elasticity  of  the  honeycomb.  In  the  calculations, 
using  program  VARYB,  the  spring  constemts  are  assumed  to  vary  linearly 
with  strain  (strain**rate  dependence  was  neglected)  and  values  were 
selected  to  provide  the  appropriate  time  duration  of  the  fall  pulse. 

In  figure  11a,  shot  1703,  the  calculated  U1bU2»31  fps  and  the 
corresponding  crush  time  obtained  fr<XB  the  experimental  U1(T)  is  780 
usee.  Good  agreement  of  A1(T)  with  e)qperiment  for  the  required  crush 
time  is  obtained  for  (F0,VI)  *=  (80,000,  0.5). 

To  indicate  the  extent  to  whi^»  the  dynamic  c2rush  force  varies 
with  the  rate*of-8train,  the  requited  crush  time  is  shown  in  figure  11a 
for  (100,000  / 0).  In  this  case,  the  calculated  Al  are  much  smaller  than 
those  of  the  es^erixssnt  of  the  beginning  and  too  large  toward  the  end  of 
the  emsh.  In  particular,  the  crush  area  is  constant  and  the 
experiments)  dAl/dt<0  for  T>440.  The  calculated  Al  can  only  decrease 
with  decreasing  dynamic  crush  force. 
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with  reference  to  the  change  in  crush  tine  with  FO  and  vi# 
calculations  show  that  *dFO/dt>constant  ■ +170  Ib/psec  in  the  range 
75,000^F0<?5,000  lb  and  -dVl/dT:0.0045/usec  for  0.4<Vl<0.5.  For 
precision  in  determining  the  crush  time  to  within  10  '»isec.  FO  is 
determined  to  within  2000  lb  and  VI  Is  determined  to  within  0.05. 

Stress  wave  effects  giving  rise  to  etbrupt  changes  in  A1  are 
apparent  at  tines  80«  420/  and  46C  psec«  However/  since  the  velocity  at 
the  end  of  the  crush  is  independent  of  F,  apd  since  the  e3q>erimental 
Ul(T)  iS/  in  effect/  a time-wise  integration  of  kl,  the  crush  time  is 
little  affected  by  stress  waves  of  short  duration*  The  calculated  data 
do  not  include  stress  wave  effects;  otherwise  the  A1(T)  agreement  with 
e^^rinent  is  within  about  10  percent* 

The  e^^riment  shows  A1>0  at  T«8S0  (+0  and  ->25)  ysec.  The 
calculations  give  Al**0  at  834  and  856  usee  for  C1bC2«0.004  and  0.008 
in.«  respectively. 

Similarly/  experimental  and  calculated  data  are  shown  in 
figure  lib  for  shot  1709.  Tlie  calculated  U1»U2*36  ft/sec  and  the 
corresponding  crush  tine  obtained  f r^  the  es^erimental  U1  (T)  is 
630  usee.  The  required  crush  time  is  obtained  for  (FO/V1)«(82/000/0.$). 
Here,  dFO/dTlconstant  « 220  Ib/usec  for  Vl»0»5  in  the  range 
80,000<FQ5.90,000  lb. 


The  e;q>eriment  shows  that  A1»0  at  1*700  ysec  -and  then  A1  abruptly 
rises  and  falls  to  zero.  The  abrupt  rise*  is  probably  ccused  by  a stress 
wave,  and  will  be  discounted.  The  fall  pulse  wets  calcuUted  for  the  two 
cases/  C1*C2*0.008  and  0.012;  the  resulting  total  pulse  times  are  692 
and  706  utec. 

5.2  Conical  Shaped  Hitigator  on  rhe  Bird  with  the  Spring  MEM 

In  the  case  of  the  spring  MEM,  the  end  of  the  mitigator  crush 
cccurs  at  the  instant  Ul*l^.  Define  UC  an  the  velocity  marking  the  end 
of  the  crush.  For  given  values  of  the  system  masses,  UO,  and  if 
UI*U2*j6sUC,  Qcmentum  conservation  yields  a unioue  value  for  UC  that  is 
independent  of.  F iind  F3.  Then,  as  in  the  previous  section,  the  time 
duration  of  the  mitigator  crush,  TC,  is  determined  from  the  experimental 
Ul^T).  With  this  value  for  TC,  we  can  proceed  to  select  a set  of 
(FQ,V1)  pairs  and  corresponding  calculated  data  to  fit  the  ei^erimental 
A1(T! . 

However^  for  the  spring  KEM,  depending  on  (F0,V1)  and  F3,  it  is 
possible  for  U2|U6  at  the  instant  Ul*U6.  Accordingly,  UC  is  not 
wiquely  de+ermined  from  momenttim  conservation.  Furthermore,  because  of 
elasticity  in  thn  honeyconb  and  KEM,  there  is  no  marked  reduccion  of  A1 
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at  T»TC,  and  consequently#  TC  cannot  be  precisely  determined  from  the 
eaqperimental  A1(T).  Nevertheless#  for  TC  obtained  from  the  condition 
U1*U2«U6»UC#  calculations  showed  that  small  variations  of  (F0#V1)  had 
little  effect  on  Al  amd  Ul  for  all  T<TC#  whether  or  not  U2bU6  at 
Consequently#  in  our  procedure#  TC  was  found  for  a (F0#V1)  pair  yielding 
U1bU2bU6«*UC;  then#  for  the  same  TC#  (FO#Vl)  pairs  were  selected  for  use 
in  the  calculations  so  as  to  fit  the  e3q>erimental  Al  for  0<T£TC.  All 
calculations  were  made  using  the  programs  PULSBl  and  PULSE3  (sects.  4.2 
and  4.3),  which  assume  linear  and  nonlinear  honeycomb  elasticity# 
respectively.  As  was  the  case  for  the  mass  MEM,  the  beginning  of  the 
crush  of  the  mitigator  can  usually  be  determined  to  within  30  ysec. 

A series  of  tests  for  evaluating  fuze  components  were  run  in  the 
HDL  4-ln.  gun,  all  with  U0s460  ± 2,  Ml»1480#  M3“260,  the  spring  MEM  of 
section  2#  and  the  conical  mitigator  shown  in  figure  4.  Here#  the 
measured  static  crush  was  107,000  t 10,000  lb.  The  crush  area#  A# 
varied  linearly  with  the  bird  penetration  (Yl-^YO)  in  the  manner 

A • 12.56  (Yl-Y6)/2.85  for  0^Y1*Y6£2.85  and 

A » 12.56  for  Yl-Y6>2.85 


For  the  given  test  conditions#  the  calculated  bird  penetration  (Y1-Y6) 
was  less  than  2.85  in.  The  measurement  error  in  the  determination  of  Al 
is  random  and  amounts  to  40Q  g for  film  readings  of  the  bird 
displacement  at  50  vsec  intervals.  The  measurement  error  changes 
approximately  as  the  inverse  square  of  the  film  reading  tine  interval. 
The  measurement  error  is  400  g for  the  given  data#  except  that  the  error 
IS  200  g for  the  data  given  in  figure  16b, 

The  esq^erimcntal  Al(T)  for  the  seven  tests  are  shown  in 
figures  12a  and  13a  where  each  test  is  identified  by  the  prefixed  letter 
E.  The  scatter  of  the  experimental  Al  is  generally  within  ±10^  g.  Some 
of  the  scatter  arises  from  stress  waves#  which  can  be  seen  from  the  data 
as  sudden  changes  in  oAl/dT.  Some  scatter  also  arises  from  measurement 
error  and  the  error  in  deteneining  the  beginning  of  the  crush.  For  scxne 
mitigatoi  shapes  with  large  altitudes#  and  also  depending  on  the 
honeycc«ib  crush  strength#  the  crush  does  not  con?>letely  occur  in  a 
column**like  manner.  In  this  case#  the  resulting  dynamic  crush  strength 
is  somevdiat  less  theui  expected  for  the  given  cross-*sectional  area. 


Momentum  conservation  gave  Ul*U2*U6«UO*85  ft/sec  and#  following 
the  above  procedure#  the  Ul(T)  from  the  seven  tests  gave 
T0665  i 10  usee.  This  suggests  that  the  error  in  determining  the 
beginning  of  the  crush  for  che  ed>ove  tests  was  closer  to  ±10  usee.  In 
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Figure  12a,  Oosparison  of  experimental  and  calculated  setback 
acceleration  of  7 shots  for  conical  mitigator 
attached  to  bird  with  spring  MEM.  Shots  1850) 
1852-1857,  Ml-1480,  M3-260,  U0»460,  and  T0665. 

addition,  it  is  interesting  to  observe  that  the  e)q>erimental  velocity  at 
the  end  of  the  pulse  for  the  seven  tests  was  -6  ± 2 ft/sec.  Apparently, 
the  e:q>eriments  are  repeatable  With  a high  degree  of  precision; 
generally,  A1  is  repeatable  to  within  tS  percent. 

The  calculated  A1(T)  for  TC  of  665  and  683  and  various  (F0,V1)  are 
given  in  figures  12b  and  13b,  respectively.  For  each  TC,  the  agreement 
between  the  calculations  of  each  set  is  ±1000  g,  except  toward  the  end 
of  the  pulse  where  the  spread  Increases  to  5000  g.  Conqparing  the  two 
sets  of  calculations  at  VI  of  0.7,  0.5,  and  0,  the  maximum  difference  in 
A1  at  any  T<665  amounts  to  1000  g.  The  calculations  for  VI  of  0.3  and 
0.2  are  also  shown  in  figures  12a  and  13a,  respectively.  Excluding  the 
fall  part  of  the  pulse,  agreement  with  the  experimentaL  data  is  about 
±1000  g.  For  T0665,  good  fit  with  the  experimental  data  is  obtained 
for  the  (FO,  VI)  range  (117,000,  0,2)  to  (102,000,  0.5).  Once  again,  as 
can  be  seen  from  the  curve  for  Vl-0,  the  dynamic  crush  force  is  a 
function  of  (U1-U6)  and  is  significantly  larger  than  the  static  crush 
force. 
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Figure  I2b.  Calculated  setback  acceleration  for  conical  mltigator 
attached  to  bird  with  spring  HEM  for  various  dynamic 
crush  strengths.  TCs»665. 

Assuming  linear  honeycomb  elasticity,  program  PUZ/SEl  calculations 
for  the  fall  part  of  the  pulse  are  shown  in  figures  12a  and  13a  for 
Cl«0.06  eind  C2«0«01,  where  Cl  and  C2  are  the  displacements  of  the 
honeycomb  at  the  bird  and  KEH  ends,  respectively,  arising  from  relaxing 
the  loads  thereon  occurring  at  T«TC*  Hence,  the  selection  of  Cl  and  C2 
'and  the  forces  acting  at  each  end  of  the  honeycomb  determine  the  spring 
constants  for  the  honeycomb;  their  values  were  chosen  so  that  the 
calculations  would  provide 

(1)  the  experimental  terminal  Ul,  and 

(2)  the  experimental  pulse  fall  time. 

Calculations  showed  that  the  experimental  terminal 
Ul*-6  ± 2 ft/sec  would  be  obtained  for  some  values  of  Cl,  C2  when  Cl  + 
C2  B 0,065  ± 0.01.  However,  in  order  to  avoid  large  oscillations  of 
Al,  which  do  not  appear  in  this  experimental  data,  calculations  show 
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Figure  I3a.  Comparison  of  experimentaX  and  calculated  setback 
acceleration  of  7 shots  for  conical  mitigator 
attached  to  bird  with  spring  MEM,  Shots  1850; 
1852-1857,  M1»1480,  H3=260r  U0=460,  and  T0688. 

that  C1»C2«  For  the  above  Cl  and  C2  limits,  the  calculated  A1  largely 
exceeded  those  of  the  experiments  in  the  neighborhood  of  T>TC. 

Consequently,  the  matching  of  the  experimental  Ul  resulted  in  calculated 
pulse  fall  times  that  were  always  too  small  (figs,  12a  and  13a). 

Apparently,  the  model  of  three  linear  springs  is  not  correct,  and 
allowance  for  elasticity  involving  the  equivalent  of  some  nonlinearity 
in  the  springs  is  required  to  provide  an  adequate  pulse  fall.  In 
particular,  the  A1  response  is  very  sensitive  to  the  honeycomb 
elasticity. 

Table  lA  is  a summary  of  a typical  calculation  of  the  forces 
acting  on  and  the  motion  of  the  various  cort^nents  of  the  system  for 
(F0,V1)  ^ (121,000,  0,3)  and  (122,000,  0,3)  for  several  values  of 

C1,C2,*  The  accelerations  Al,  A2,  A6  are  given  in  units  of  10^  g and 

the  forces  F,R,F3  are  in  units  of  10^  lb.  The  bird  penetration  of  the 


*ln  the  calculations,  TC  » 665  at  (F0,V1)  « (121,150,  0.3). 
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Figure  13b.  Calculated  setback  acceleration  for  conical  mitlgator 
attached  to  bird  with  spring  MEM  for  various  dynanac 
crush  strengths.  TO=688. 

honeycomb,  (Y1-Y6),  amounts  to  2.0  in.,  at  which  point  the  crush  force 
is  88,000  lb.  For  each  FO,  the  maximum  spring  MEM  loading  of  86,400  lb 
IS  reached  between  500  and  600  usee.  It  is  interesting  to  note  that  the 
hydrodynamic  force  R»M4(U1-U6)  arising  from  the  honeycomb  crush  amounts 
to  15  percent  of  F at  lOG  usee.  Thereafter,  R/F  falls  rapidly,  since  R 
varies  as  (Ul-U6)2.  Although  TC  differs  by  20  ysec  for  the  two  values 
of  FO,  the  change  of  A1(T)  is  negligible.  Of  course,  because  of  the 
same  UC,  the  difference  in  A1(T)  between  the  two  pulse  fall?  is 
acco'vinced  for  by  this  20  ysec  difference.  For  each  Cl,  C2,  we  observe  a 
rapid  fall  of  A1(T)  with  oscillation. 


In  PULSES,  Z1«Z2*-A1*M1/X1  and  A3'*Zl* (Xl-X2)/M3,  where  A1  is  the 
experimental  A1(T)  and  XI,  X2  are  the  calculated  Insteuntaneous  honeycomb 
elongations  arising  r..om  the  instantaneous  force  at  each  end  of  the 
honeycomb  at  T^TC.  In  this  way,  the  honeycoasb  elasticity  is  made 
nonlinear.  The  spring  MEM  elasticity  is  linear  and  is  the  same  as  in 
PULSEl.  In  PULSE3,  the  displacements  Xl«X2*C3»C2  at  T«TC  are  selected 
input  to  provide  tne  same  pulse  fall  time  as  t .at  of  the  experiment. 


Using  PULSE3/  the  calculated  results  for  (FO,Vl)s(i21,000/  0.3}  emd 
several  values  of  C1«C2  are  shown  in  table  IB.  Clearly,  the  time 
required  for  X3  to  go  to  zero  is  essentially  independent  of  the  selected 
Cl,  C2«  Moreover,  X3-0,  at  the  same  time  the  experimental  Alto.  Hence, 
the  spring  MEM  determines  the  time  duration  of  the  pulse  fall.  Finally, 
C1bC2  are  uniquely  determined  by  the  condition  that  simultaneously 
X1-X2-X3“A1**0  at  the  end  of  the  fall. 

For  the  values  shown  in  table  IB,  the  fall  time  increases  with 
increasing  Cl»C2.  For  the  total  pulse  time  of  1.02  msec,  to  within 
10  percent,  C1»C2*0.08  is  the  required  honeycomb  elongation  at 
T»T0674  msec.  Excluding  the  neighborhood  Al>0  where  the  procedure  is 
doubtful,  Z1bZ2  is  largest  near  T°TC  and  decreases  with  increasing  T by 
a factor  of  about  3. 

Excluding  the  fall  part  of  the  pulse,  the  present  calculations 
agree  with  those  of  the  previous  section  for  the  mass  MEM,  where  for  the 
measured  static  crush  of  90,000  ± 9000  lb,  a good  fit  was  found  for 
(F0,V1)  * (80,000,  0.5).  However,  the  shape  of  the  present  mitigator 
differs  from  that  used  for  the  mass  MEM,  and  this  can  result  in  a change 
of  VI  in  characterizing  the  dynamic  crush  strength. 

5.3  Wedce  Miticator  on  the  Bird  with  a Spring  MEM 

Two  tests  were  run  in  the  HDL  4-in.  gun  at  U0«453  with  Ml*1630, 
M3*230,  the  spring  KEM  of  section  2,  and  the  double  wedge  mitigator 
shewn  in  figure  5.  The  altitude  and  base  lengths  were  l.S  and  0.9  in., 
respectively.  The  crush  area  varied  approximately  linearly  with  the 
bird  penetration  (V1-Y6),  so  that 


A » 12.56(Y1-Y6)/1.5  for  Yl-Y6<1.5  and 
A » 12.56  for  Yl-Y6^1.5  . 

Momentum  conservation  gave  UO90  and  the  TC  from  the  experimental 
Ul(T)  were  found  as  550  and  580  ysec  for  shot  numbers  1897  and  1898. 
The  e3q>erimental  terminal  velocities  differ  by  only  2 ft/sec.  Again, 
li)ce  those  for  the  conical  shaped  mitigator,  the  present  experiments 
appear  to  be  very  repeatable  and  the  beginning  of  the  crush  can  be 
determined  to  within  30  ysec. 

The  calculated  Al(T)  for  TC-564  with  (FO.Vl)  « (86,000,  0.5)  is 
shown  in  figure  14  for  T<TC.  The  scatter  of  the  experimental  data  as 
well  as  the  fit  with  the  calculations  is  generally  within  ±1000  g.  The 
details  of  the  forces  acting  on  and  the  motions  of  the  various 
coiqponents  are  summarized  in  table  II.  The  bird  penetration  (Y1-Y6) 
slightly  exceeds  the  l.S-in.  altitude  length  of  the  wedge. 
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Figure  14.  Comparison  of  experimental  and  calculated  setback 

acceleration  for  double  wedge  mitigator  attached  to 
bird  with  spring  HEM.  Ml=1630,  M3«230,  U0«=453,  and 
TC«564.  Shots  1897  and  1898. 

The  pulse  falls  of  the  two  tests  are  in  close  agreement,  but  no 
attempt  was  made  to  calculate  the  pulse  falls  or  determine  the  spring 
constants.  The  slightly  higher  experimental  terminal  velocity  found  for 
tests  1897  and  1898,  5±1  ft/sec,  con^ared  with  those  of  the  previous 
section  (-612  ft/sec),  is  due  in  part  to  the  slightly  higher  efficiency 
of  the  gun  at  higher  values  of  (M1+M3),  (which  results  in  higher 
pro3ectile  energy  and  a higher  value  for  UC)  and  may  also  be  due  to  a 
very  small  difference  of  honeycomb  elasticity  between  tlie  two 
natigators. 

Similar  agreement  (not  presented  here)  in  experimental  scatter  and 
between  experimental  and  calculated  data  is  found  for  the  triple-wedge. 
The  Al(T)  for  the  two  types  of  wedges  are  also  in  good  agreement. 

5,4  Mitigator  Placed  Adjacent  to  Spring  MEM 

It  is  preferable  to  place  the  mitigator  adjacent  to  the  MEM, 
rather  than  accelerate  it  in  a gun,  since  the  entire  gun  energy  can  then 
be  used  to  drive  the  bird  and  payload.  For  a given  maximum  setback 
deceleration,  this  leads  to  increase  of  UO  and  a longer  pulse  time. 
Again,  to  avoid  strong  stress  waves  and  obtain  a smooth  pulse,  the 
mitigator  crush  must  begin  at  the  shaped  face.  For  this  reason  the 
shaped  face  is  turned  around  so  that  now  the  crush  proceeds  from  the 
bird  face  toward  the  MEM.  As  noted  in  section  3,  this  affects  the 
forces  acting  on  Ml  and  M6. 


The  double-wedge  mitigator  described  in  the  previous  section, 
figure  5,  with  an  altitude  of  1.5  in.  and  a base  length  of  0.9  in.,  was 
placed  with  its  base  adjacent  to  the  spring  MEM.  In  this  test,  HDL  shot 
nuitiber  1901,  U0*464,  Ml»1800,  and  M3*210. 


Momentum  conservation  gave  UC»88  ft/sec  and  TC  from  the 
experimental  U1(T)  was  found  to  be  €00  usee.  Again,  the  time  of  the 
beginning  of  the  crush  is  known  only  to  within  about  30  usee,  so  that 
calculations  were  made  for  TC  of  598  and  632  usee  at  (F0,V1)  of 
(82,000,  0.5)  and  (77,000,  0. 5) , respectively.  To  form  a conraon  basis  of 
comparison,  3C  psec  was  subtracted  from  the  plot  of  the  latter 
calculations.  As  shown  in  f. gure  15,  the  agreement  between  the 
calculated  emd  experimental  data  is  generally  to  within  2000  g for 
T0598  usee  and  about  1000  g for  TC=632  usee.  The  details  of  the  forces 
acting  on  and  the  motions  of  the  various  coR5>onents  are  summarized  in 
table  III  for  TC*632. 


« 00  ro  VI  uc  TC  iioeo.  $ 4m 
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Figure  15.  comparison  of  experimental  and  calculated  setback  accel- 
eration for  double-wedge  roitigator  attached  to  spring 
MEM.  Ml*1800,  M3=210,  U0=4b4,  TC=59B,  and  TC*632. 

Shot  1901. 
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5.5  Additxone.1  Data 

Additional  tests  were  performed  with  double-wedge  mitigators  for 
static  crush  strengths  of  725,  2X30,  and  4350  psi  at  bird  impact  speeds 
between  309  and  650  ft/sec.  The  results  of  these  tests  and  calculations 
bassd  on  the  conservation  equations  are  given  in  figures  16,  17,  and  18. 
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Figure  16a  shows  7-in.  gun  test  data  and  the  results  of 
calculations  using  VARYB  for  00*648,  P*725,  and  with  the  mitigator 
ittached  to  a mass  MEM.  Test  data  were  not  recorded  for  T>0.7  msec,  so 
that  part  of  the  pulse  during  the  mitigator  CAish  and  the  entire  pulse 
fall  is  missing.  Figure  16b  gives  4-in«  gun  test  data  and  the  results 
of  calculations  using  VARYB  for  U0*318,  F*725,  and  with  the  mitigator 
again  attached  to  a mass  MEM-  For  each  test,  the  difference  of  Al(T) 
between  experiment  and  calculation  for  (F0,V1)  * (1.2P,0.3)  is  well 
within  10  percent  of  A1(T).  The  fall  pulse  shown  in  figure  16b  arises 
entirely  from  the  elasticity  in  the  uncrushed  honeycomb.  As  noted  in 
section  5.3  and  in  the  present  calculation,  a good  fit  with  the  test 
data  is  obtained  for  C1*C2=0.01.  Here,  the  pulse  fall  amounts  to  about 
350  Msec,  compared  with  the  100  ysec  fall  obtained  in  section  5,1.  The 
reason  for  the  longer  fall  arises  from  the  larger  values  of  Ml  and  M2 
and  smaller  A1  at  the  time  U1*U2«UC. 

The  shot  of  figure  16t  was  designed  to  obtain  a nearly  constant 
pulse  of  long  duration  by  controlling  the  mitigator  area.  Accordingly, 
the  double-wedge  altitude  amounted  to  2 in.,  and  the  mitigator 
cross-sectional  area  linearly  increased  in  this  length  from  0 to  only 
7,8  in.^  This  area  then  further  increased  linearly  to  12.56  in.^  in  a 
lenqth  of  16  in.  The  experimental  bird  penetration,  6,  amounted  to 
17,5  in.  and  the  predicted  value  was  17.8  in.  Generally,  the  difference 
between  calculated  and  test  values  of  6 for  either  bird  or  MEM 
penetration  of  the  mitigator  is  within  eUx>ut  4 percent  of  6. 

Calculated  and  test  data  A1(T)  for  450<u0<500  with  P of  2130  and 
4350  and  the  mitigator  attached  to  the  MEM  arc  shown  in  figures  17  and 
18,  respectively,  A better  fit  with  the  test  data,  especially  the  slope 
dAl/dT  from  the  maximum  value  of  A1  to  A1  at  T=TC  is  given  by 
(F0,V1)  « (1.05P,  0.5)  than  for  (1.2P,  0,3).  Again,  for  each  set  of 
calculations,  the  difference  between  the  calculated  and  test  data  Al(T) 
IS  within  10  percent  of  Al(T), 


Apparently,  the  dynamic  force  (eq  12)  adequately  describes  the 
crushing  of  honeycomb  of  various  strengths  in  the  range  725;^P£8000  psi 
at  crush  speeds  up  to  at  least  650  ft/sec.  In  agreement  with  the 
physical  observation  of  the  crushes,  this  means  that  honeycomb  crush 
occurs  in  approximately  the  same  m«mnQr  over  the  above  ranges  of  P and 
crush  speed. 
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Figure  16a*  Comparison  of  experimental  and  calculated  setback  accel- 
eration in  the  7-in*  gun  for  double-wedge  honeycomb 
nitigator  attached  to  mass  MEM*  Wedge  altitude  « 5 in., 
P=725,  Ml=2530,  M2=31280,  M3=2160,  UO=64B,  and  (FO»Vl)  ^ 
(1.2  P,  0.3).  Shot  53.  (Data  not  available  for  T>0.7  msec.) 


Figure  16b.  Comparison  of  experimental  and  calculated  setback  accel- 
eration for  double-wedge  honeycomb  mitigator  attached 
to  infinite  mass  HEM.  Wedge  altitude  « 2 in. , 

P»725*  Ml“4180,  M3«1220,  UO^SIS,  Cl»C2*=0.01  and 
(FO,Vl)«(1.2  F,,  0.3;.  Shot  1984. 
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Figure  17,  Con5>arison  of  experimental  and  calculated  setback  accel- 
eration for  double-wedge  honeycOTib  mitigator  attached 
to  bird  with  mass  MEM.  Wedge  altitude  * 1,5  m,, 

P«2130,  Ml«1580,  M2"30430,  M3«320#  U0*454,  and 
C1»C2»0,01,  Shot  2022, 

6.  SUMMARY  AND  CONCLUSIONS 

The  use  of  variously  shaped  and  various  strength  aluminum 
honeycomb  mitigators  together  with  stock-item  railroad  springs  in  the 
MEM  has  demonstrated  that  it  is  possible  to  obtain  smoothr  predicted>le, 
highly  repeatable  setback  acceleration-time  pulses  for  testing 
projectiles  and/or  components  contained  therein.  The  experimental 
scatter  from  test-to-test  is  found  to  be  1000  g for  accelerations 
attaining  a maximum  of  25,000  g«  The  difference  between  experimental 
and  theoretical  setback  accelerations  based  on  the  mass,  momentum,  and 
energy  conservation  equations  is  less  than  10  percent  for  the  rise  and 
steady  parts  of  the  pulse.  The  crush  can  usually  be  predicted  to  within 
4 percent  of  the  bird  or  MEM  penetration  of  the  mitigator.  The  fall 
pulse  is  governed  by  the  nonlinear  system  elasticity  and  the  physical 
constants  thereof  arc  determined  from  experimental  data. 
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’Figure  18.  Coaparison  of  e^g^rimental  ai^  calculated  setback  accel* 

eration  for  double-wedge  honeycomb  aitigator  attached  to 
bird  with  spring  MEM.  Vtedge  altitude  * l.S  in.# 

P-4350,  Ml=1250,  m3«300,U0=480,  C1»0.06,  and  C2»0.0l. 
Shots  2014  and  2015. 

The  dynamic  force  equation  (12)  is  derived  from  the  experimental 
data  and  is  of  the  same  form  for  mitigator  strengths  of  725  to  8000  psi 
and  crush  speeds  up  to  at  least  (>50  ft/sec.  The  difference  between 

experimental  and  calculated  setback  accelerations  is  within  10  percent 
for  either  (FO, VI)  pairs  (1.05P,  0.5)  or  (1.2P,  0.3).  Pulse  shapes  may  be 
varied  over  a wide  range  of  rise  and  steady  times,  peak  acceleration, 
etc. , by  varying  the  mitigator  shape  (crush  area  as  a function  of  crush 
length),  crush  strength,  projectile  mass  and  spring  MEM  elasticity. 

The  spring  MEM  stores  energy  during  t))e  crush  of  the  mitigator 
and,  together  with  some  small  energy  provided  by  elasticity  of  the 
mitigator,  releases  this  energy  at  the  termination  of  the  crush.  For  an 
energy  corresponding  to  a linear  spring  loading  of  86,400  lb  at  a 
deflection  of  0.331  in.,  the  spring  MEM  provides  a smooth  pulse  fall  of 
about  300  psec  for  a 3.8  lb  projectile  attaining  a peak  setback 
deceleration  of  25,000  g.  This  compares  to  a sharp  pulse  fall  of  about 
100  psec  for  a mass  MEM,  where  the  elasticity  is  derived  solely  from  the 
uncrushed  mitigator. 
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7.  LIS®  0?  SYKBOLS 


A mltig&tor  crash  erea  at  the  bird  or  MEM  interface,  in,^ 

Ai  acceleration,  ft/scc^ 

Cl  honeycomb  elongation  at  the  bird  interface,  arising  from 

relaxing  the  force  thereon  at  T«TC,  in. 

C2  honeycomb  elongation  at  the  HEM  interface,  arising  from 

relaxing  the  force  thereon  at  *i>TC,*in. 

C3  displacement  of  the  spring  in  the  spring  MEM  at  T«TC,in. 

EO  either  the  initial  kinetic  energy  of  the  bird  (mitigator 
attached  to  HEM) , or  initial  kinetic  energy  of  the  bird  plus 
edtigator  (mitigator  attached  to  bird) , ft- lb 

£1  energy  dissipated  by  mitigator  deformation,  ft- lb 

E2  instantaneous  kinetic  energy  of  the  system,  ft-lb 
T mitigator  dynamic  crush  force,  Ib 

FO  product  of  dynamic  crush  coefficient  and  P,  used  to  determine 
mitigator  dynamic  crush  force,  equation  (12) , lb 


F3  force  exerted  by  the  spring  in  the*  spring  MEM,  lb 
I»1  bird 


1»2  mass  hEM  or  aft  section  of  spring  MEM 
1=3  original  uncrushed  mitigator 


I«4  crushed  mitigator 

I«S  remaining  uncrushed  mitigator 

1=6  forward  section  of  spring  KEM 

t distance  of  mitigator  crush  front  travel,  in. 
(t-6)  length  occupied  by  crushed  mitigator,  in. 


L original  mitigator  length,  in. 
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list  of  symbols  (Cont'd) 

HO  mass  of  railroad  springi  g 

MI  masst  g 

H4  (»pAS{Ul-U6)  for  spring  HEM;  =pAS (01-'J2>  for  mass  MEM), 
rate  of  nitigator  crush,  iWsec 

p nitigator  static  crush  pressure,  psi 

R (.  K4(Ui-«6)  for  spri.e  MEM,  - H4(Ul-02)  for  mass  HEM), 
hydrodynamic  crush  force,  Ifc 

S {-f/S) . ratio  o^ush  front  travel  to  depth  of  bird  or  MEM 
penetration 

T time,  sec 

TC  time  duration  of  the  mitigator  crush,  sec 

UC  velocity  at  the  termination  of  the  mitigator  crush,  ft/sec 

01  velocity,  ft/sec 
00  initial  bird  velocity,  ft/sec 

VI  constant  used  in  determining  mitigator  dynamic  crush  strength, 
equation  (12),  nondimensional 

W energj'  dissipated  through  collision  between  crushed  and 
uncrushcd  nitigator  masses,  ft- lb 

YI  displacement#  in* 

« depth  of  bird  or  MEM  penetration,  in. 

p density  of  uncrushed  nitigator,  Ibn/in, 

pi  density  of  crushed  mitigator,  Ibm/in.^ 


> ,teg*Sgt*-  ■•■■rurg.-y.Cti  -t.  ^ 

. ’-__ 


LIST  OF  SYMBOLS  (Cont'd) 


For  Hass  KEH 

X1«C1-Y3+Y120,  honeycofflb  elongation  at  bird  interface,  in. 
X2=C2-Y2+Y3£0,  honeycomb  elongation  at  MEM  interface,  in. 


Z1(=-A1H1/C1) , honeycomb  spring  constant  at  bird  interface,  where  a1  is 
the  acceleration  at  T=TC,  Ib/in. 


22(*-AlHl/C2) , honeycomb  spring  constant  at  MEH  interface,  where  A1  is 
the  acceleration  at  T=TC,  Ib/in. 


For  Spring  hem 

X1=C1-Y3+Y1S0,  honeycomb  elongation  at  bird  interface,  in. 
X2»C2-Y6+Y3ao,  honeycomb  elongation  at  MEM  interface,  in. 
•C3-V2+Y6iO,  spring  displacement  of  spring  MEM,  in. 


Z1  (-A1M1/C1,  linear  case;  — AIMl/xi,  nonlinear  case),  honeycomb  spring 
constant  at  bird  Interface  where  Al  is  the  acceleration  at  T*TC» 
Ib/in. 


Z2  (=-AlHl/C2,  linear  casej  — AlMl/Xl,  nonlinear  case),  honeycomb  spring 
constant  at  MEM  interface,  where  Al  is  the  acceleration  at  T=TC, 
Ib/in. 


Z3  {=86,400/0.3308),  spring  constant  for  spring  MEH,  Ib/in. 
Superscript 


(•)  denotes  time  differentiation  of  the  giver,  variable 
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irtifi  iTiniAininintnjnino  00  ■oo'n 
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8.3  (CONT'D) 

6/0  r-T+li 

6a0  IF  U6<J!  GjrO  290 
690  GOTO  640 


' • V 


i •- 
<NO< 
OH  I O 


► C OUl  > 

- — .S£  H 

3 U f<) 

^ * O 

I 3 <i 

-j-o'n  j_  I-.  H H cj  — fvj  ' 
3r3xHa:5^2'32>*«s^  : 
i n n 9 n — «a—  m n » n 
— ttazcc  — — c\jn  > 
s'coHy.  ^a.*-^a.>-3:3=>; 


COOOOOOOOOOOOOOCOOOOOOOOC( 
c — n^>»'0oc^o<r>’«r>oor^'0^c  — 

fs.f^r«r«-p*fv.*^r»'--*>.C0‘D'»'X'DC0‘»X't)XO*O.0‘O>>C 


040  A I *- . /5*G*  ( BO-^B  I ♦Tv32*T''24B3*r'3+B4*T^44-85*T^5^B6*T*^6  7+B8*T''8 ) 

050  XI»C1-Y3^Y1 
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